Long-range magnetic fields in the ground state of the Standard Model plasma 
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In thermal equilibrium the ground state of the plasma of Standard Model particles is determined by tempera- 
ture and exactly conserved combinations of baryon and lepton numbers. We show that at non-zero values of the 
global charges a translation invariant and homogeneous state of the plasma becomes unstable and the system 
transits into a new state, containing a large-scale magnetic field. The origin of this effect is the parity-breaking 
character of weak interactions and chiral anomaly. This situation can occur in the early Universe and may play 
an important role in its subsequent evolution. 



It is generally believed that the ground state of the Standard 
Model at high temperatures is homogeneous and isotropic. 
This assumption underlies the description of all the important 
processes in the early Universe: baryogenesis, cosmological 
phase transitions, primordial nucleosynthesis, etc. fl44il In this 
work we demonstrate, however, that at finite density of lepton 
or baryon numbers due to parity-violating nature of the weak 
interactions this homogeneous "ground state" becomes unsta- 
ble by developing a long-range magnetic field. The transition 
to the "true" ground state may depend on the details of the 
non-equilibrium dynamics, when various violent dissipative 
processes (e.g. turbulence, radiation emission, finite conduc- 
tivity of plasma) play an important role. 

What are the conditions for the translational invariance to 
be spontaneously broken by a long-range field? It is suffi- 
cient for the free energy of the gauge fields to contain an in- 
teraction term that dominates over the kinetic energy and can 
be both positive and negative. An example is provided by a 
Chern-Simons term /cs oc ADA, that has less derivatives that 
the kinetic term {dAY and therefore can dominate over it at 
large scales. The presence of the Chern-Simons term in the 
Maxwell equations is known to lead to an instability and gen- 
eration of magnetic fields. 

At zero temperatures and densities the Chern-Simons term 
for electromagnetic fields is prohibited as a consequence of 
gauge invariance and Lorentz symmetry {Furry theorem 
At finite temperatures and densities the plasma creates a pre- 
ferred reference frame and the 4-dimensional Lorentz invari- 
ance is broken down to 3-dimensional one. As a result the free 
energy of static gauge fields is 

T[A] = I d^pA,{m,j{p)Aj{~p) + 0{A^) (1) 

with the polarization operator 

n,y (p) = {p'^S,, - p,p, )Ui + ze,,fcp^n2 (p2) , (2) 

where i,j,k = 1,2,3 are spacial indices; p^ = e^-fe 
is the antisymmetric tensor. Eq. ^ is the most general form 
ofHij satisfying the gauge-invariance trans versality condition 
PiUij = 0. In the long wavelength limit p^ ^ Q a non-zero 
112(0) means that the Chern-Simons term 112(0) A ■ V x A 



appears in ([T]l. The 3x3 matrix (|2|l has then a negative eigen- 
value for sufficiently small momenta p < |n2(p^)/ni(p^)| 
and the corresponding eigenmode grows larger and larger (un- 
til the higher order in A terms would stabilize it). In the above 
consideration it is important that the gauge field is Abelian. 
Unlike the Yang-Mills fields ^ the magnetic component of 
the photon field does not get screened in plasma Jst] (i.e. 
ni(0) remains finite) and therefore the instability does not 
require large 112(0). 

In this work we demonstrate that in the Standard Model 
plasma in the Higgs phase an equilibrium value of 112 (0) for 
electromagnetic fields is non-zero and proportional to the val- 
ues of the global charges: baryon (B) wd flavor lepton num- 
bers La (index a runs over flavours). Unlike the previous 
works iJ-'Sl (see discussion below) it is important that even if 
the anomalous charge B + L is absent, 112 (0) remains non- 
zero and magnetic fields develop ll45ll . 

Chern-Simons term and axial anomaly. The origin of the 
112 term has its roots in the axial anomaly (see e.g. |@, IH IH- 
[Tlll ). Indeed, the non-conservation of the axial current at finite 
densities of left or light fermions n l , nn means that one can 
convert fermions into gauge field configurations with a non- 
trivial Chern-Simons number Nqs = / (Px A ■ B (where B = 
V X A is a magnetic field): 

'±1^ = ^ fd^.E.B = ^'-^ (3) 
at 27r^ J -n at 
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(here E = —A is an electric field and a = is the 
fine-Structure constant). Let us consider the simplest exam- 
ple of left and right fermions at zero temperature with dif- 
ferent Fermi energies (chemical potentials) pL^ 7^ hr- In- 
finitesimal change of the gauge field 5 A will destroy (create) 
SriL ji = ±^SNcs of real fermions around the Fermi level. 
The total energy of the system will change by ST = (/i^ — 
f-F()^SNcs ig], which leads to the parity-odd Chern-Simons 
term in the free energy: T[A] = "^''-^^"^"^ J d^xA- B. 

This remains true in any vector-like gauge theory at finite 
temperature/density where there is a difference of chemical 
potentials of left and right-chiral charged particles [4, 9, l3l 
(also ifisi [14I1 ). Indeed, to calculate the polarization opera- 
tor we need to analyze one-loop contribution from charged 
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FIG. 1: Polarization operator [(a)] one-loop weak corrections to it ( |(b)| and [(c)] ) and its expansion in A/i/r [(d)[ 



fermions described by the diagram [Ta] If the left and right 
fermions have different chemical potentials such that 
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(where a;„ = Tr{2n + 1)T, ?? e Z are the Matsubara frequen- 
cies andPL.R = 5(1=1=75) are chiral projectors) their contribu- 
tions to this diagram are different. Assuming first that A/i <C 
T, let us consider linear in A/i/T, correction to the polariza- 
tion operator (for this one should differentiate the fermions 
Greens function © with respect to fi). This correction is de- 
scribed by the diagram[Td] with A/i playing the role of a third 
external field. The diagram [Ta] thus turns into the famous tri- 
angular graph for the axial anomaly iITsI - IitIi . with the third 
vertex containing "axial vector field" = (5^o^M7o75- The 
resulting term in the effective action, oc £aPtJ.uXaApdfj,A^, 
again reduces to the Chern-Simons term with 112 (0) = ^^M- 
This expression for 112(0) is actually exact in A/i and j[4[]. 

Similar logic applies to the non-Abelian gauge fields ^ Ul 
[ol-ITTll. In the Standard Model with its chiral weak charges of 
fermions the coefficient in front of the SU(2) Chern-Simons 
term can be expressed in terms of hb+l (B being baryon and 
L lepton numbers) QIIll. However, in this case a homoge- 
neous state becomes unstable only at large values of chemical 
potential, exceeding the mass of weak bosons. Even at high 
temperatures in the symmetric phase the "magnetic screen- 
ing" effect m requires A/i > T to overcome the "mag- 
netic mass" rriniagn ^ awT. Moreover, anomalous non- 
conservation of B + L current drives the coefficient of cor- 
responding Chern-Simons term to zero |[l9ll and the stand- 
ing wave-like configurations of the gauge fields are actually 
metastable (see discussion in [7]). 

An example of the situation when the non-zero coefficient 
in front of the Chern-Simons is realized in the Standard Model 
at high temperatures, T ^ nif, when the smallness of the 
electron's Yukawa coupling makes the number of right elec- 
trons conserved at classical level. If the initial conditions have 
non-zero /ig^, the Chern-Simons term for the U(l) hyper- 
field is then generated (with 112(0) oc /ien) the genera- 
tion of long wave-length magnetic fields occurs ([sj], until the 
chirality-flipping reactions, suppressed as (me/T)^ do not de- 
stroy the /ig^ . In the early Universe where such a situation can 
be realized, the rate of these reactions becomes comparable to 
the Hubble expansion rate at high temperatures T ^ 80 TeV. 




FIG. 2: Fermi corrections to the self-energy of the fermion. 
The loop gives non-zero contribution only at finite lepton and 
baryon number density. 



Magnetic fields, generated in such a way are rather short- 
wavelength (much smaller than the horizon size at that epoch) 
and are probably erased during the subsequent evolution due 
to plasma dissipative processes. At lower temperatures all 
the chirality-flipping reactions are in thermal equilibrium and 
naively the chemical potentials of all left and right-chiral par- 
ticles are equal. However, it was shown in |[20ll that if strong 
helical magnetic fields are initially present in the plasma, then 
the relaxation rate both for A/i for electrons and for helical 
fields significantly increases and they both can survive down 
to T ^ 10 McV. The ground state that the system eventually 
reaches contains neither fields nor A/i. 

We demonstrate below that although these considerations 
are true for electrodynamics, in the Standard Model where 
fermions are also involved in parity-violating weak interac- 
tions, the difference of chemical potentials of all left- and 
right-chiral fermions is actually present (with all chirality- 
flipping reactions taken into account) and leads to the gen- 
eration of magnetic fields. 

In this paper we analyze the simplest situation when this 
effect is present: the case T ^ mw (mass of the -boson) 
when weak interactions can be described by the Fermi theory; 
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The full Hamiltonian of the theory H = Hq + Hp + H em has 
a free part for fermions and photons. Ho, and terms describing 
electromagnetic ("Hem) and Fermi ("Hf) interactions |[46ll . 

Dispersion relation of fermions and chemical potentials. 

To describe the equilibrium plasma at T ^ mw we intro- 
duced the density matrix, g ~ Z^^ cxpl — /?('H — XaLa — 
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^qQ^^bB) j . Five global charges commute with the Hamil- 

tonian H: three La, B and Q (Aq, Aq, As are the correspond- 
ing Lagrange multipliers); and the partition function Z en- 
sures that tr(^) = 1. 

To find the distribution functions of left and right-chiral par- 
ticles we compute the correlators (?/'Pl,r'/') = 
We expand the density matrix in interactions to get 



(6) 



At zeroth order in interactions one gets {t/jPLt/j 



('0Pr'0)o = 5('0'0)o- The propagators of left and right par- 
ticles (up to corrections of the order m/T ^ 1) have the 
form dUi with fi^ — M-R- This conclusion remains true if we 
take into account parity-preserving electromagnetic interac- 
tions. 



Taking into account chiral Fermi interactions, gp w ^o(l ~ 
PHf) one finds that ('0PlV')f (V'^r'0)f- The Green's func- 
tion is found via = Gq^ — S where for example the 



self-energy of left electron, Egt , is 



4Gf 



(7) 



Expression (O does not depend on momentum, the thermal 
averages {ip'tp)o and {iP^^Pl.r'iP)o are proportional to the 
particle-antiparticle asymmetry (see e.g. MM)- To compute 
(5/iei one should substitute in ^ the thermal averages, sum- 
marized in the Table 1 in Appendix B. As a result, for example 
the electron propagator becomes 



G = 



1 



(8) 



i.e. the dispersion relation of electrons change when taking 
into account Fermi corrections (c.f. ll2l1l . ll47l '). Indeed, from 
(|i — TO — = we see that the "on-shell conditions": 
(w — /itree)^ ~ + gcts shifted for left (right) particles 
by 2{uj — pixee)5iiL,R (in the limit 5ijll,r ^ ^^), where ^tree = 
(Aq — Ae). In the Hmit rrie/T — >■ 0. Eq. ([8]) splits into the 
sum of free propagators in the form (0) with different chemical 
potentials Pl.r- 

This difference can give rise to a parity-odd term in polar- 
ization operator of photons 101. Indeed the polarization oper- 
ator that was parity-even when computed with respect to the 
density matrix Qq acquires a parity-odd part when averaged 
with respect to the gp. The lowest order weak corrections 
are represented by two diagrams, [Tc]and [lb] The computation 
of the diagram [lb] is quite similar to that of |4] and gives a 
non-zero n2(0) = ^ - S^ij^): 



n2(0) 



a AGp 



2^ y/2 



CL^La + CbB 



(9) 



where coefficients cl^,cb ^ C'(l) depend on the fermionic 
content of the plasma (details are summarized in Appendix) 
ll48ll . The diagram [Tc] does not contribute to the 112(0) as it 
can be cut into two diagrams[Ta]along the vertical dotted line, 
each of which is at least first order in momentum (Eq. ([2)). 

Although the Fermi theory ^ is not renormalizable, the 
result ([9]l is given by the non-divergent part of the diagram [2] 



and is expressed in terms of well-defined physical quantities 

(c.f. iH). 

Chern-Simons coefficient at two loops and "non- 
renormalization tlieorems" . The diagram [id] is similar 
to the triangular diagram, responsible e.g. for 7r° — > 27 
decay (with A/i?/'7o75V' playing the role of the only non-zero 
component of the chiral current, describing pion) iHsI - IitI . 
It is well known that the axial anomaly should be calculated 
at one loop only and is not renormalized by higher-loop 
corrections ll22l - l25ll . also at finite temperature and density. 
At the same time our result becomes non-zero only at 
two loops. There is, however, no contradiction. What is 
non-renormalized for the chiral anomaly is the numerical 
coefficient in front of the proper combination of external 
fields, (e.g. ^ in Eq. In our case this coefficient is also 
not renormalized. The structure of the parity-odd one loop 
term has the same form at tree-level and at one-loop in G^: 
112(0) = ^(A/i,ree + (5/i)A • B. The numerical coefficient is 
dictated by the axial anomaly; AfitKe is a possible difference 
of chemical potentials present at tree-level (zero in our case); 
and 6p is the shift generated by the diagram[Tb] The structure 
of the 112(0) term therefore remains the same as in [0] with 
the total difference of chemical potential (A/iu-ee + Sp). 

Also a four dimensional theory at finite temperatures can be 
regarded as a three dimensional Eucledian model albeit with 
the infinite number of particles - each Matsubara mode of a 
fermion becomes a "particle" with the mass a;„ = 7r(2n+l)r, 
n G Z. Therefore (as it was argued in um) our result 
may seem to be in contradiction with the "Coleman and Hill 
(CH) theorem" 1I27I1 that states that in any Eucledian three- 
dimensional gauge theory without massless particles 112 (0) = 

l^J^^i exact at one loop. However, the pres- 

ence of the infinite number of modes changes the situation, 
as can be seen already in the simplest chiral gauge theory, 
if one computes the 112(0) in the Matsubara formalism (see 
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e.g. 101). Formally, considering the left- and right-chiral par- 
ticles as fermions with "complex mass" to„ = (cj„ — i/iL.j?)- 
and applying directly to results of ll27ll one would arrive 
to an undefined expression 112 (0) = §- Yl 

Uri-ifj-R rpj^g reason why this happens is clear: the 

degree of divergence of the diagrams is different in 3 and 
4 dimensions (hence the infinite sum over n). In particu- 
lar, if we first summed over the Matsubara frequencies and 
then integrating over momentum (or if one uses dimensional 
regularization in 3-momentum integral and then takes a limit 
d ^ 3, iHl, see also ifTlll ) one obtains a well-defined answer 

of H 

Moreover, the CH theorem uses the fact that the 3-point 
vertex r^^\pi . . . ) — 0{pi). This is not true in our case, 
as the diagram [id] becomes linearly divergent in 4 dimen- 
sions and therefore the shift of integration momentum by any 
fixed vector k changes its parity-odd part by a finite amount 
(X Afie'^^-AjknAj Smtl. 

Ground state. The presence of 112(0) ^ leads to the 
generation of magnetic fields. The Chem-Simons number 
A'cs ~ kA^ will increase until it reaches f-A^cs ^ {nL — 
nil) ^ GpLtot (see e.g. i^). At fixed A^cs the magnetic field 
tends to increase its wavelength to decrease the total energy 
(B^ ~ fciVcs). As a result, the system does not have a ther- 
modynamic (infinite volume) limit (c.f. j6;]): the value of the 
field and the scale of the inhomogeneity will be determined 
by the size of the system. It is clear, however, that in realistic 
systems establishing of the long-range field is a complicated 
process (see e.g. i20ll ). greatly affected by the dissipative pro- 
cesses and by existence of different relaxation channels of iVcs 
(resistivity of plasma, energy radiation, turbulence, etc., see 
e.g. ll30l - D2ll ). This may significantly affect subsequent evolu- 
tion and even the final state of the system. 

Discussion. In this work we demonstrated that the Standard 
Model plasma at finite densities of lepton and baryon numbers 
becomes unstable and tends to develop large scale magnetic 
fields. We considered electrodynamics plus Fermi theory (|5), 
a description of weak interactions valid when e^™" < 
(T/mw)^, i.e. at T < 40 GeV. At higher temperatures one 
should consider full electroweak theory and perform two-loop 
computations of 112(0). At even higher temperatures (in the 
symmetric phase) one should analyze hypermagnetic fields. 
We leave these analyses for future works. We expect however 
that our conclusion about the instability of a homogeneous 
state will hold. 

Below we discuss several realistic systems in which the ef- 
fects discussed here can become important. As a first exam- 
ple, let us consider the primordial plasma at radiation domi- 
nated epoch. Eq. (|9]l gives 112(0) ~ c x a{GFT^)r]L,B where 
Vl.b < 1 is the ratio of the total lepton (baryon) number to the 
number of photons, = ^''^P and the numerical coeffi- 
cient c «a 2.5 X 10"^. We see that the CS coefficient decreases 



with temperature fast and therefore the effect is the strongest 
at high temperatures T < mw ll49ll . 

The instability starts to develop at scales k ^ 112(0) and 
the magnetic field initially growth as where p = k^tja 
(see e.g. IM [20(1). The conductivity of the plasma is ct ^ 
O(10^)T 133D. The requirement for an instability to develop 
over the characteristic time of temperature change (i.e. Hub- 

ble time) is: /3(T) « 2.0 (;;^) ' (f§4)' > 1. We see 
that the measured value of baryon asymmetry ijb ~ 6.0 x 
10^^" |[34ll is too small to trigger any instability. The sit- 
uation is different for lepton asymmetry where only the up- 
per bounds at the epoch of primordial nucleosynthesis exist: 
^ few X 10~^ 1I35I1. At earlier epochs even 77^ 1 is 
possible (if this lepton asymmetr y di sappears later). This is 
the case e.g. in the z^MSM (see Qql for review), where the 
lepton asymmetry keeps being generated below the sphaleron 
freeze-out temperature and as a result may reach the levels 
rjL ^ 10^2^ 10"^ before it disappears at T - few GeV iH. 
We see that significant magnetic fields can develop in this 
case, which can play an important role for analysis of the cos- 
mological implications of the i^MSM. 

As a next application we consider a high density degenerate 
electron plasma (appearing e.g. in white dwarfs and neutron 
stars, OSll ). Notice, that our consideration remains valid in 
this regime, as Eq. ([7]i makes no assumption about the rela- 
tion between mass, temperature and chemical potential of the 
particles. Only the numerical coefficient in Eq. ^ changes 
and we checked that it is non-zero and 0(1). The same re- 
lation 112(0) ^ j^GpLtot holds, however now Ltot = n^ 
(density of electrons) that can be quite essential, reaching 
10^" 10^^ cm"3 in the crust of neutron stars iHl. The 
corresponding scale of the instability k ^ 112(0) is then in 
(sub)km size and the time of its development is much shorter 
than the lifetime of the star. 

To summarize: in this work we discussed a previously un- 
known effect that occurs in the Standard Model at finite tem- 
perature and density. It implies that a number of processes 
in the early Universe can be affected, including cosmolog- 
ical phase transition, baryogenesis, dark matter production. 
This effect may in particular lead to the generation of horizon- 
scale helical cosmic magnetic fields purely within the Stan- 
dard Model. Such fields may survive till present and serve as 
seeds for the observed magnetic fields in galaxies and clusters. 
The effect may also be important for explanation of physics of 
compact stars. 
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Charged leptons: 








{U) = i7o(-An,) 


Neutrinos: 








{^^) - \ 


PL7o(-Anj.) 


Quarks: 


{q^Pq) 


= <5'^°iAn5 






70 (-An,) 



TABLE I: Thermal averages of fermions. Particle-antiparticle asymmetry is defined as An = f / ^J'yi [/p — /p] , fpifp) is the 
Fermi-Dirac distribution for particles (anti-particles) and the number of internal degrees of freedom, .g = 2 for neutrinos v; 

5 = 4 for charged leptons i; g ~ 12 for quarks q. 



Appendix A: Fermi theory 



For completeness we summarize in this Appendix the definitions of charged and neutral currents in the Fermi theory (see e.g. 
Chapter 20 in ll40ll '). The charged currents are defined as 

J^'^'+ = {e-Lltj^VL + dLlt^UL) + other generations , (Al) 
v2 

and neutral currents are 

i' 

Here Pl,r = ^(liTs) ^re chiral projectors, charges = {T^ — sin^ dwQ), g% ~ {— sin^ dwQ), T3 = ±i is the 3rd generator 
of SU{2), Q is the electric charge and the sum in (IA2b goes over fermions in all flavours. 



Appendix B: General expression for the Chem-Simons coefficient 112(0) 



A general expression for 112 (0) is given through the asymmetries of all fermions (if some fermions are absent in the plasma, 
their asymmetry should be put to zero). 



n2(o) = 



a AGf 



2tt ^/2 
1 

^ 324 
67 



I E " § ~ 2 cos{29w)) J2 



9 ^ ° 36 

a 

17 - 62 cos{29w)) {^riu + Auc) 
(91 + 134 cos(26'vy )) {And + Aus) 



(Bl) 



+ ^(l + 2cos(26lM/))Anfc 

(here 9w is the Weinberg's angle). 

Once the asymmetries of all particles An are expressed through the conserved charges B and La under the condition of 
electric neutrality of the plasma ll43ll , the expression (IB It reduces to the form (|9): 



n2(0) 



27T y/2 



CL^La +CbB 



(B2) 



where the values of coefficients depend on the fermionic content of the plasma. For example, if plasma contains 5 quarks (except 
for the top quark) and all leptons, then 



C-Lo 



22cos(26'iv)-45 
621 



CB 



53cos(26'w) +430 
621 



(B3) 



